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TITLE: AUTOMATED OPTICAL MEASUREMENT APPARATUS 
AND METHOD 



CROSS REFERENCE TO RELATED APPLICATION 

This application is a Continuation-In-Part of United States Patent Application 
No. 09/328,972 filed on June 9, 1999 and also claims priority from United States 
Provisional Patent Application No. 60/160,949 filed on October 22, 1999, both of 
which are incorporated by reference herein in their entirety. 

BACKGROUND OF THE INVENTION 

This invention generally relates to methods and apparatus for testing optical 
systems and particularly to devices and techniques for the automated measurement of a 
variety of parameters of optical surfaces and/or elements including radii of curvature, 
surface shape, thickness, power, focal length, wavefront, topography, and aberrations. 

Throughout the process for fabricating optical systems or components, it is 
frequently necessary to determine if, and how well, optical surfaces or elements 
conform to their designers stated requirements. Not only does the performance of 
optical systems in final form need to be verified, but various parameters of their 
components need to undergo intermediate testing for conformance with their 
specifications prior to final assembly as a system. Indeed, even the tools of fabrication, 
especially molds for the formation of plastic or glass lens elements, need to be tested 
for compliance with design specifications. Some of the most frequently encountered 
measurements that need to be made are radius of curvature of surfaces in either convex 
or concave form, surface topography, thickness, power, and various focal lengths. 

Classically, radius of curvature is measured through the use of a hand-held 
instrument called a spherometer, which measures the sagittal height (sag) of the surface 
over a known diameter and then displays the radius of curvature on a dial or other 
visual display after an internal calculation that relates radius to sag height and the 
known diameter. However, the accuracy of such devices are prone to relatively large 
errors because sag heights are usually small dimensions that are difficult to accurately 
measure mechanically. 
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A more accurate technique for radii measurement involves the use of an auto- 
collimating microscope in an arrangement referred to as a radiusscope. Here, one first 
focuses on the surface to be measured and then on the center of curvature of the 
surface where a reticle image has been imaged back on itself by reflection from the test 
5 surface. The positions of the microscope are recorded, and the difference between 
them represents the radius of curvature to limits of accuracy, which depend on the 
preciseness of the length measurements and the ability of the operator to accurately 
focus. 

Where the spherometer suffers from problems of precision, the use of the 

10 radiusscope, which can be accurate to microns if care is taken, is time consuming and 
dependent on operator skill and experience. 

The thickness of an optical element is more or less important depending on its 
assigned role in a particular design and can be critical where the design relies heavily 
on its precision for aberration control or the like. Thickness obviously can be 

15 measured directly by mechanical means, which may also be automated, but there is 
always the danger of damaging part surfaces with mechanical approaches. 

Power and focal length are always of interest and can be calculated from 
classical lens makers formulae having knowledge of the various numerical values 
required as, for example, index of refraction, radii, and thickness. 

20 For rapid qualitative tests other approaches are also often employed. For this 

purpose, reliance on the Foucault knife-edge test is often made, and while the Foucault 
test has its uses, it suffers from an inability to sense changes in slope deviation that are 
changing only slightly in magnitude or direction. For more accurate quantitative 
results, reliance is made on interferometric and screen methods, such as the Hartmann 

25 test, to perform wavefront analysis from which more gross properties may be derived. 
However, such tests, while informative, have been laborious and time consuming. 

While those skilled in the art have provided a variety of ways for measuring 
many of the foregoing properties of optical elements and systems, there remains a need 
for an instrument for rapidly and accurately measuring a number of optical properties 
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virtually simultaneously, and it is a primary object of the present invention to provide 
such a device. 

Another object of the present invention is to provide methods and associated 
devices for automatically performing wavefront analysis, measuring surface 
5 topography, and determining radius of curvature, thickness, power and focal length of 
optical surfaces and/or elements with minimal dependence on operator skill. 

Yet another object of the present invention is to provide an automated 
instrument for providing statistical analysis of quality in high volume production 
settings. 

10 Another object of the present invention is to provide apparatus and methods for 

ri optically measuring the physical characteristics of surfaces that are at least partially 

reflective. 

% -\ Still another object of the present invention is to satisfy all of the foregoing 

Q1 objects with a user-friendly device that is simple in its implementation and low in cost. 

,ll 15 Other objects of the invention will in part be obvious and will in part appear 

f hereinafter. A full understanding of the invention will best be had from reading the 

Q detailed description to follow in connection with the detailed drawings. 
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SUMMARY OF THE INVENTION 



The invention comprises a system and method for automatically performing 
dynamic screen testing on a surface and determining its shape from which other optical 
parameters of interest may be derived and reported. A measuring head, consisting of a 
source, beamsplitter, objective lens, and lens array with a CCD camera, is mounted on 
25 a translation stage and moves along the optic axis of the head relative to the part under 
test. The part under test is mounted on an appropriate support, such as a three-point 
support nest, that automatically centers spherical parts on the optical axis of the 
system. 

Light is projected along the optical axis through a microscope objective or 
30 other appropriate lens to illuminate the part under test with a predetermined wavefront, 
preferably spherical, so that subsequent calculations are made simpler if this light is 
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recollimated parallel to the optical axis of the system. Light reflected from the part 
under test passes back through the lens, after which it passes through a pellicle or cube 
beamsplitter towards a CCD camera. A two-dimensional array, preferably in the form 
of a pair of crossed lenticular screens, is placed in front of the CCD active area so that 
5 a series of sharp images are formed on the CCD array. When the system measuring 
head is positioned so that the focal point of the objective is located near the surface of 
the part under test or near its center of curvature, then the incoming nearly parallel 
light produces a series of spots on the CCD active area. The shifts in the pattern of 
spots are used to determine the shape of the surface under test. Mathematical analysis 

10 of this shape provides information on the radius of curvature of the part (if spherical), 
the "Spherical" and "Cylindrical" radii of curvature of a toric part (along with the 
angle between the major axes and a given reference axis), and the "Shape Factor" of an 
aspheric part. For ease of interpretation, the overall shape can be expressed in various 
ways, including Zernike polynomials. Software performs this analysis and facilitates 

15 providing results in many useful forms - contour plots, wire- frame models of 

deviation, direct readout of coefficients, direct readout of RMS surface form, direct 
readout of peak-to-valley difference, etc. Display screens are customizable for the 
engineering specialist or on-the-floor auditing and measurement for production. 
Custom processing capabilities are available using Visual Basic® and an Object 

20 Linking And Embedding (OLE®) interface. 

DESCRIPTION OF THE DRAWINGS 

The structure and operation of the invention together with other objects and 
25 advantages thereof may best be understood by reading the detailed description in 
connection with the drawings wherein parts have been assigned a unique reference 
numeral that is used to identify them whenever they appear in a drawing and wherein: 
Fig. 1 is a diagrammatic perspective of the automated measuring apparatus of 
the invention showing an enclosed optical head in combination with a computer 
30 workstation; 
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Fig. 2 is a diagrammatic front perspective view of the optical head of Fig. 1, 
uncovered so that its major components may be seen; 

Fig. 3 is a diagrammatic, right quarter perspective view of the optical head of 
Fig. 1, uncovered so that its major components may be seen from a different vantage 
point than that of Fig. 2; 

Fig. 4 is a diagrammatic, left quarter perspective view of the optical head of 
Fig. 1 uncovered so that its major components may be seen from a vantage point 
different than that of either Fig. 2 or Fig. 3; 

Fig. 5 is a diagrammatic elevational view of the major components of the 
optical head of Fig. 1 showing how a spherical test surface is illuminated with a 
predetermined wavefront originating at the surface's center of curvature and how that 
wavefront is reflected from the test surface and subsequently measured; 

Fig. 6 is a diagrammatic elevational view of the major components of the 
optical head of Fig. 1 showing how a spherical test surface is illuminated with a 
predetermined wavefront focused on the surface and how that wavefront is reflected 
from the test surface and subsequently measured; 

Fig. 7 is a diagrammatic elevational view of the major components of the 
optical head of Fig. 1 showing how a spherical test surface is illuminated with a 
predetermined wavefront originating at some location other than the surfaced center of 
curvature and how that wavefront is reflected from the test surface and subsequently 
measured; 

Fig. 8 is a diagrammatic illustration showing the anticipated location of an 
array of dots from a spherical surface illuminated with a spherical wavefront 
originating at the surface's center of curvature along with deviations of the dots as the 
origin of the wavefront passes through the surface's center of curvature (through 
focus); 

Fig. 9 is an enlarged diagrammatic perspective showing how "crossed 
lenticular" screens are used in the invention to provide an array of equivalent spherical 
lenses; 




Fig. 10 is a screen from a program used in conjunction with the apparatus of the 
invention to provide a graphical user interface by which properties of a test surface can 
be displayed, various system instructions issued, and other test unit properties 
described; 

5 Figs. 1 1 A-D are output screens from a computer program used in conjunction 

with the apparatus to graphically display topographic maps of the surface of a test unit 
for various states of defocus; 

Fig. 12 A-D are output screens from a computer program used in conjunction 
with the apparatus of the invention showing various surface contour plots for different 
1 0 states of defocus of a test unit; 
P Figs. 13 and 14 are output screens illustrating various coefficient plots as a 

function of the position of a wavefront origin with respect to a test surface for a toric 
surface and spherical surface, respectively; 
pi Fig. 1 5 is an output data screen showing measured values for a toric surface; 

15 Fig. 16 is a diagrammatic elevational view similar to Fig. 5 with the addition of 

3 a lens forward of the microscope objective to enable measurement of long radius parts; 

H Fig. 17 is a diagrammatic elevational view similar to Fig. 5 except for the 

absence of the objective which has been removed to enable measurement of flat 
P surfaces that are at least partially reflective; 

20 Fig. 18 is a diagrammatic elevational view similar to Fig. 17 except for the 

presence of an additional mirror and a relay system for measuring the wavefront error 
of transmissive filters over their bandpass; 

Fig. 19 is a diagrammatic elevational view similar to Fig. 17 with the addition 
of a beam expansion section to measure large flat surfaces; 
25 Fig. 20 is a diagrammatic elevational view similar to Fig. 18 except for the 

absence of the relay system and measures the wavefront error of transmissive filters 
over their bandpass; and 

Fig. 21 is a diagrammatic elevational view of an alternative embodiment of the 
invention that uses an intervening telescope between the objective and CCD array to 
30 make full use of the available numerical aperture (NA) of the objective. 
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DETAILED DESCRIPTION OF THE INVENTION 

The invention comprises a system and method for automatically performing 
dynamic screen testing on surfaces of optical components or systems or other surfaces 
that are at least partially specularly reflecting. The inventive system and method may 
5 be used for rapidly and accurately measuring optical parameters of a variety of optical 
elements and/or surfaces including, but not limited to, those of lenses and mold 
surfaces for forming them, flats, filters, and prisms. Lenses may be spherical, aspheric, 
or torics. The system is capable of measuring wavefronts, surface topography, radii of 
curvature, thickness, power and focal length, wavefront error in bandpass filters, and 
10 can exist in more or less complex forms depending on a selection of available features 
j=i and desired range of operation, but in all of its forms, its principle of operation is the 

~l same. 

y - 

Within less than a minute, for example, a measurement of radius of curvature 

m can be made to an accuracy of less than 2.0 micrometer error, and other parametric 

15 measurements can be carried out with similar speeds. One form of the inventive 

s system with the versatility of automatically measuring concave and convex elements 

p over a wide parametric range and along different azimuths of rotation is shown in Figs. 

y 1 where it is designated at 10. As can be seen in Fig. 1, along with Figs. 2-4, automatic 

i y 

Q measuring system 10 comprises a vertically oriented optical head 12, a Pentium based 

20 computer 14, a monitor 16, preferably color, and a keyboard 18 that is preferably an 
enhanced type. Optical head 12 is electrically connected to computer 14 via an 
interface board (not shown) located in one of computer 14's expansion slots. The 
interface board is configured in a well-known manner as an I/O board to: provide an 
electronic link between computer 14 and optical head 12 for the transmission of 
25 electrical and logic signals between them; provide drive signals to stepper motors 
located in optical head 12; and perform certain signal processing operations on raw 
signals generated by a video camera in optical head 12. Miscellaneous housekeeping 
chores may also be implemented via the arrangement. A separate framegrabber board 
may also be resident in one of the expansion slots of computer 14. Suitable cables (not 
30 shown) provide electrical links. 
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Keyboard 1 8 operates in the usual way to provide the operator with a means for 
communicating commands to computer 14, and monitor 16 displays a variety of 
system responses to the operator including measurement data, system status, and menu 
options to be selected for system control and direction. This is preferably done via a 
5 graphical user interface (GUI) in conjunction with a mouse 20 or other suitable device. 

In addition, data and computer instructions in the form of programs may be 
transferred between operator and computer via one or more disk drives such as those 
shown typically at 34 in Fig. 1 . The inventive system may also be networked in well- 
known manners. 

10 On top of optical head 12 there is preferably provided a rotary stage 38 for 

receiving and supporting test optical elements in a reference plane so that they may be 
yP measured. 

•vj A fixed reference plane RP, as best seen in Fig. 5, is defined by three nylon 

^ balls 40 (only two of which are shown) that are precisely located in a V-groove formed 

15 in a nest ring 42 which also has a tapered conical tunnel that extends through to the 
3 interior of optical head 12 and is centered about its optical axis, OA, although not 

~l shown. In this connection, it is to be noted that the three nylon balls 40 also force the 

Q center of curvature of spherical surfaces to lie along the optical axis, OA, of the 

p instrument. This is so because the various optical components comprising the optical 

~" 20 head 12 are precisely aligned with the optical axis which, in turn, is made parallel to 
the mechanical axis of the optical head 12, and the mechanical axis is preferably made 
precisely perpendicular to the plane defined by the balls 40. If only spherical surfaces 
are to be measured, it is also possible to use a donut instead of the three nylon balls 40. 
As can be seen in Fig. 5, test elements, such as that shown as 26, 
25 advantageously are automatically located in the reference plane by this three-ball 

location arrangement, even if slightly tilted about their own axis. When viewed from 
above, a circle centered on optical axis, OA, and drawn through the centers of each of 
them defines the semidiameter of the reference plane. 

Referring now to Figs. 2-4, which show the interior of optical head 12, it is 
30 seen to further comprise a scanning head 50, which can be vertically moved along 
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optical axis, OA, via a translation stage 56 that, in turn, is driven via a conventional 
lead screw (not shown). The precision lead screw is driven by a precision stepper 
motor 58 fixed to the rear wall of optical head 12 and under command of computer 14. 
Translation stage 56 moves parallel to optical axis, OA, on a complementary 
5 configured precision slide arrangement not shown in detail, but generally designated at 
62 in Fig. 4. The pitch of the lead screw and the rotational steps of motor 58 are 
preferably selected in a well-known manner so that each step of the stepper motor 58 
controllably moves the scanning head 50 up or down by a small fraction of a 
millimeter. The accuracy of this stepping motion can be enhanced by the use of 

10 encoders which directly measure translation of stage 56 at each step. In addition to 
improving accuracy, encoders can also be utilized in minimizing inaccuracies due to 
wear, especially in the precision lead screw assembly 62. 

Computer 14 is programmed in a well-known manner to keep track of the 
position of scanning head 50. In addition, limit switches may be provided to prevent 

15 scanning head 50 from mechanically interfering with any structure beneath rotary table 
38 and also to provide a locating signal indicating that scanning head 50 has reached 
its uppermost position when closed. In this manner, means are provided for precisely 
moving scanning head 50 along optical axis, OA, while at the same time providing a 
signal by which its vertical position can be monitored and controlled via computer 14. 

20 Scanning head 50 comprises a number of components all of which travel 

together along optical axis, OA. As best seen in Fig. 5, with occasional reference back 
to Figs. 2-4, these comprise a microscope objective 70, an adjustable diaphragm 71, 
preferably a pin hole, a collimated source of illumination 74 preferably in the form of a 
light emitting diode (LED), provided with suitable collimating optics, and a pellicle 

25 beamsplitter 80 located between objective 70, source 74, and a video camera 76 to 
provide a split optical path between the three. The objective for camera 76 is 
preferably a pair of crossed lenticular screens with index mismatching material 77 
sandwiched between them as best seen in Fig. 9. The crossed lenticular screens 
operate to provide the system with an array that is equivalent to a 2D array of spherical 

30 lenses each of substantially the same focal length. The presence of the index material 



-9- 



provides a mechanism by which the focal length of the individual equivalent spherical 
lenses of the screen may be varied by a factor of four or more depending on the 
availability of the index of the material in accordance with the following relationship: 



where «' is the index of refraction of the lenticular screen material and n is the index 
of refraction of the material 77 between the crossed lenticular screens. Here, common 
indices range from approximately from 1.40 to 1.59. 

It will be recognized that the lenticular screen may be replaced by an array of 
spaced apart mechanical apertures or holes, but with an attendant reduction in the 
dynamic range of the apparatus. This is so because the spherical lenses make for more 
tightly focused spots thereby making the amount of change in the local slope greater 
before spots overlap compared with what would be the case with pure mechanical 
apertures. Also, use may be made of commercially available microlens arrays in place 
of the crossed lenticular screen and intermediate mismatching material. 

Filters may also be used as needed to control the spectral output of the LED. 
Light source 74 preferably comprises an LED although a laser, fiber optic source, 
incandescent lamp, or discharge lamp may be used. Considerations in the choice of the 
type of source to use include efficiency, cost, lifetime, ease of adjustment of light 
levels, and compactness. The spectral output of the source 74 should be considered 
because objective lens 70, being a microscope objective, is in the best case corrected 
for visible wavelengths. So, the spectral content of source 74 needs to be considered in 
terms of its impact on signal level because of any longitudinal aberrations introduced 
at wavelengths outside the visible region. In addition, it is important to note that the 
part being measured may introduce chromatic aberrations if measured with a broader 
bandwidth source, and this may influence measurements in an undesirable way. 
Consequently, an LED is quite acceptable because it is a narrow source with no 
speckle or laser coherence problems. In addition, it is fairly monochromatic and easy 
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to modulate. The bandwidth of a narrow LED source for this purpose is typically tens 
of nanometers wide, not hundreds, and such devices are readily commercially available 
at low cost. In addition, they require low power to operate and have low heat 
dissipation. 

5 The magnification, aperture, and NA of objective 70 are selected in well- 

known manners to operate in conjunction with the other system components so that 
adequate signal level is available over the desired area of the surface to be measured. 

Having described the major system components of optical head 12, its 
operation will now be. taken up referring to Figs. 5 -7. Optical head 12, comprising 
10 source 74, beamsplitter 80, objective lens 70, and lens array 78 with CCD camera 76, 
is mounted on translation stage 56 and moves along the optical axis of the head relative 
to the part under test. The part under test is mounted on an appropriate support, such 
as a three-point support nest that automatically centers spherical parts on the optical 
axis of the system. 

15 Light in the form of a collimated beam with a plane wavefront 75 is projected 

along the optical axis via beamsplitter 80 through microscope objective 70 or other 
appropriate lens to illuminate the part under test with a predetermined wavefront, 
preferably spherical. The spherical wavefront is preferred, although not strictly 
required, so that subsequent calculations to be described are made simpler if the 

20 reflected light is re-collimated parallel to the optical axis of the system. Light reflected 
from the part under test passes back through microscope lens 70, after which it passes 
through pellicle beamsplitter 80 towards array 78 of CCD camera 76. Lenticular array 
78 is placed in front of the CCD active area, one focal length away from it, so that a 
series of sharp images are formed on the CCD array. This image or series of images as 

25 the case may be, is passed via framegrabber 77 to computer 14 for analysis. 

When the system measuring head is positioned so that the focal point of the 
objective is located near the surface of the part under test as shown in Fig. 6, or near 
the center of curvature of the part under test as shown in Fig. 5, then the incoming 
nearly parallel light produces a series of spots on the CCD active area , Ax = 0 and 

30 Ay = 0, the x andy displacements measured as deviations of a spot's centroid (See Fig. 
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8 which shows the variation of spot locations with focus position). The shifts in the 
pattern of spots is used to determine the topography of the surface under test as 
described hereinafter. 

As best seen in Fig. 7, when the origin of the spherical wavefront is at some 
5 point other than the surface or its center of curvature (assuming a spherical surface), 
the return wavefront is no longer a plane wavefront and Ax ^ 0 and Ay * 0 . As will be 
seen , mathematical analysis of these deviations leads to a measurement of the surface 
topography which, in turn, provides information on the radius of curvature of the part 
(if spherical ), the "Spherical" and "Cylindrical" radii of curvature of a toric part (along 

10 with the angle between the major axes and a given reference axis), and the "Shape 

Factor" of an aspheric part, as well as other parameters which may be readily derived. 
For ease of interpretation, the overall shape can be expressed in various ways, 
including Zernike polynomials. Software performs this analysis and facilitates 
providing the results in many useful forms as, for example, contour plots, wire-frame 

15 models of deviation, direct readout of coefficients, direct readout of RMS surface 
form, direct readout of peak-to-valley difference, etc. Display screens are 
customizable for the engineering specialist or on-the-floor auditing and measurement 
for production. Custom processing capabilities are available using Visual Basic and an 
Object Linking And Embedding (OLE) interface. 

20 Having described the components of the optical head and how plane and 

distorted wavefronts are generated in it, the mathematical algorithms used to analyze 
the information so generated will now be taken up. It is assumed that the part to be 
measured can be described, at least approximately, as having a radius of curvature r. 
Z-bar, the surface shape is then the sum of two components, a z-component related to 

25 the radius of curvature r and the sag h(x,y), which is the departure of the surface from 
the base radius, r. Thus the surface may be represented by a sum of coefficients and 
functions, preferably as algebraic polynomials. Further, the focal length of the 
objective (or other lens) used to focus the collimated beam is f. The focal length of the 
array of lenslets is F. The separation between the array and the back principal plane of 

30 the objective is d, and the separation between the array of lenslets and the CCD is 
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nominally F, but it may be different, in which case the actual separation is to be used in 
calculations in place of F. Moreover, if a known perfect spherical part is measured, the 
actual distance between the CCD and the array of lenslets may be deduced from the 
resultant data and thereafter used in subsequent measurements. 

Whenever the position of the optical head 12 is located so that the focus of the 
objective is very close to the surface of the part under test or its center of curvature , 
the reflected light will be very nearly collimated, and the result will be a series of well- 
focused spots on the CCD, with very nearly equal spacing between adjacent spots 
(assuming the lenslets in the array are nearly equally spaced). There will, in general, 
be a shift between the location of a given spot and the spot one obtains when truly 
collimated light is incident on the CCD camera with its lenslet array. So long as the 
light is nearly collimated, there is no great difficulty in associating any given deviated 
spot with the position of an undeviated spot. If necessary, one such spot can be 
decreased in intensity or otherwise identified by eccentricity in shape, intensity, or 
location so as to provide a reference point to help identify spot correspondence. (The 
appearance of a series of undeviated spots can easily be generated by placing a plane 
first- surface mirror atop the lens nest, or use a reflection from the unused "arm" of 
light reflected from the beamsplitter). The size and direction of the deviation of this 
spot from the expected position provides information on the angular shift of the ray of 
light at that point. This can be related to the difference in the part under test between 
the expected slope for a spherical surface and the actual slope. From this is derived the 
deviation in surface height between a spherical reference surface (whose radius is 
equal to the distance between the "Surface" location and the present position of the 
translation stage) and the true surface. Note that the actual surface Z-bar will always 
be substantially the same regardless of the location of the objective focus with respect 
to the surface provided sensible deviations are being measured. 

The difference in height h(Xp i j,y P ij ) of the surface of the part under test is 
expressed as a polynomial in two variables, x pi j and y P iji 

h ( x Pv>ypy) =a io x pij + a oxy P ij + a nX P ijy Pi j+ a 2o x2 P ij + a o2y Pij + ■•• 
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Here the subscript "p" denotes a coordinate on the surface of the part (without 
the "p" the coordinates x and y refer to coordinates in the CCD plane). The subscripts i 
and j refer to the numbering of the coordinates/spots. The angular deviation 8 xij at a 
spot (x p ,y p ) measured along the x-direction is given by: 



dh 



dx r 



^ xij 



Similarly, the angular deviation s yi j along the y-direction is given by: 



dh 



= £ 



{xptj>ypii) 



yy 



Straightforward calculation (using, as an example, polynomials only to second order) 
then gives: 



£ xiJ = a w +2a 20 x pij +a n y f 



The relationship between these slopes and the observed shifts of the spots seen 
by the CCD can be related by the following expressions, obtained from geometrical 
optics. The value of the radius of curvature, r, should be set equal to the separation 
between the positions of the translation stage at which the spot shifts are minimized : 
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f 

£ xij =AX vJ]^ = a ™ + 2a 20 X piJ + a ny P & 



f 

Zyij = AVy — = <*oi + a u x P ij + 2fl ra ^. 



Finally, the coordinates on the surface of the (concave) part are related to those 
in the plane of the CCD by the expressions: 



dr_ 
f 



"(r + /) 



y on 



v/y 



r 



dr r 



Here Axy and Ayy are the shifts in the x- and y- directions of the ijth spot. (For convex 
parts we replace r by -r). The index i runs from 0 to N, while the index j runs from 0 
to M. 

Next, the values of the coefficients of the polynomials are calculated by 
performing a least-square fit of the observed shifts of the spots to the shape of the 
surface, using the above relations. If e X jj = £y and e y ij = r|ij, then this least-square fit 
can be accomplished by solving the matrix equation: 




pys py 



0 

2 H x Pijy P ij 
o 



4 



*20 



V a 027 



Y&P&j +x Pi/ny) 
2 T x Pi£v 
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The relevant expressions for carrying the surface profile out to fourth order are 
longer and more daunting. The slope equations become: 

%y =a io +2a 20 x + a u y + 3a 30 x 2 + 2a 2X xy + a n y 2 +4a 40 x 3 + 3a 3l x 2 y + 2a 22 xy 2 + a X3 y l 
7] /y = a ox + a u x + 2a 02 y + a 2l x 2 + 2a X2 xy + 3a 03 y 2 + a 3X x 3 + 2a 22 x 2 ,y + 3a n xy 2 + 4a M y* 

5 

(To save space, the subscripts on the variables have been suppressed. All x's are x P ij 
and all y's are y pi j, all variables on the surface of the part.). The matrix corresponding 
to the second-order matrix above is a 14 X 14 matrix. It will not fit on this page 
legibly, so it is broken down by a vertical fissure into two halves, which are reproduced 
10 below, the left side first (Note that the matrix is symmetrical about the diagonal. 
Again, the subscripts have been suppressed for the purposes of clarity.): 
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NM 
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2 2> 


5> 


0 


3 2> 2 
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NM 


0 
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2> 2 


pi 
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Z(- 2 ^ 2 ) 


2^xy 
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Z(2*w) 


4EV 
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0 
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3^>y 






0 
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12^ 5 












z( 3 *v + * 4 ) 


2l^V 




^(6xV + x 5 ) 
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6ZV 
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The right half of the matrix is : 



" o 4jy 3j?y 2530? 2y 0 

ijxy s^y 0 jy 2^y y 35^ 4^y 

^2*W) 45y> 5f*V+* 4 ) 2 lP + ^) H^+W) 4£x>? 

4j>y 6jy 0 2^ 4^yy sjy 

3W 0 i2jy 9j>V s^yy sjyy 0 
25p+^ 35V/ 8^>v _26*y+^) xl^y+^v) jfv+^y) 4jyy 
Xf^y+y) 6^>y 4jyy ^y+^v) ^v+w) x^+e^y) s^y 4 
sjy 0 3£*y ejyy 9j> 4 i2jy 

4jyy 0 i£> 6 12^ s2>v 45vy 0 
5j3x?y+2cv) 32yy i©^ _X^y+* 6 ) 3fr?y+26) 3jW+x 4 y) 4jyy 
jjvwyj 6^yy 8jxV 26x 3 y+2x 5 >) 4y(^y+xv) iJ^+^j s^yy 
jjywy) 9j> 4 4^y 35^y+*y) ^y + ^y) jjy-f^y) 1253^ 

sj>y i2^y 0 4jyy s^yy 15}^ i£> 6 



5 The accompanying "Results" column vector is a 14 X 1 expression: 



X( 3 *>i^ + *i^) 
2 E(v^ 

Z(y^- +3 w>«) 
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As will be discussed in more detail hereinafter, the calculations can be 
extended for surfaces of any power, but it should be noted that the number of terms 
grows very rapidly. For a surface expressed in powers up to the nth order there will be 
5 (n + l)(n + 2)/2 terms, and the required matrix expression will be Q x Q, where 

Q = {(n + l)(n + 2)/2}-l. 
Once the coefficients have been derived, it is a simple matter to calculate the 
surface departure from a sphere, and to calculate the rms surface departure, the peak- 
to-valley difference, and any other measure desired. The size of the reference sphere 
10 against which the surface is compared, it should be noted, has a radius equal to the 
p distance between the "Surface" location (position at which the focus of the Objective 

lens coincides with the surface of the Part Under Test) and the present location of the 
^ focus of the objective (in general, there will be no result unless this location is either 

01 near the "Surface" location itself, or is near the Center of Curvature of the part). 

.ji 15 The departure from sphericity of the surface can be expressed in the form of 

= Zernike polynomials (which are simply appropriately weighted combinations of 

p surface polynomials), which are commonly used descriptors of optical elements, useful 

m in the study of aberrations. 

The radius of curvature of a spherical surface can be measured in several ways. 
20 The simplest and most direct is to measure the distance between the two "crossing 
points" of the second order coefficients. This distance is equal to the radius of 
curvature of the spherical surface. Provided the part is aligned so that the optical axis 
runs through the center of the spherical surface, the plots of ao2 and a 2 o as a function of 
position z will be identical, falling to zero at the same positions near the "Surface" and 
25 the "Center of Curvature" (CoC) locations. This is a very direct and robust 

measurement, since it will yield the radius of curvature even if the instrument is not 
properly calibrated. That is, the relationship between the spot shifts Ax^ and Ayy and 
the deviation of the surface from a true sphere might not be correctly calculated 
because the distances or focal lengths might not be properly entered, but (as long as the 
30 focus scale of z measurement is properly calibrated) the separation between the 
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crossing points will give the radius of curvature in all cases. In fact, this method can 
be used to check the calibration of the translation stage. 

The radius can also be obtained from the variation of the coefficients a 2 o and 
a 02 with position, z. Again, assuming that both curves coincide (so that we may speak 
5 of them as the single function, a 2 (z)), and assuming z is measured from the Surface 
position, the relationship between a 2 (z) and z will be given by: 



11/ \ 1 o 

- = - + l«20 +^0 2 J = - + 2a 2 



r z 



10 



Which can be re-arranged to give 

1 1 

a 2 = 

A plot of a 2 versus position can yield the radius in one of two ways. If a 2 is plotted 
against z, then one can derive the radius from the slope: 



15 



20 



r = 



da 2 
Ik 



A plot of a 2 vs. 1/z gives the result more directly, since the intercept, by the 
above equation, is l/2r. Both of these methods can be used to calibrate the instrument, 
if the measurement is performed with a known spherical standard. (The radius of the 
part need not even be known, since the coefficient of 1/z must be -1/2, All one needs 
to know is that the part is spherical.) 

Finally, if the part is known to be spherical, but the plots of a 2 o and a 02 do not 
overlap, then this separation can be used to calculate the decentering of the part (i.e. - 
the separation between the center as defined by the fixture for the part and the optical 
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axis of the measuring head). One can express the coefficients a 2 o, a 0 2, and an as a 
function of the center location (x 0 ,yo), or vice-versa: 



a 20 = 2 



1 1 ^ ip +y 2 Q 



l 1,1 x 2 o+yl 



a„rs 



^( a 20 - a ^f~ a \\ + (^20 " ^02 ) 



10 When the surface has a well-defined geometrical form which is nearly 

spherical, the coefficients can be analyzed to yield information about that shape. A 
good example is a toric surface, in which the surface resembles a slice from the outer 
edge of a torus. Such a surface has one radius of curvature along one axis and a 
different radius of curvature along the axis perpendicular to the first. A toric shape is 

1 5 commonly used in contact lenses to correct for astigmatism. A sphere is a special case 
of a toric in which the two radii of curvature are equal. 
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The general equation for a toric can be given in the form: 



x 2 + y 2 + z 2 +R 2 -r 2 -2R^x 2 + y 2 =0 

These variables x, y, and z are different from those used in the rest of this paper, and 
5 should not be confused with them. This equation describes a circle of radius r which is 
rotated about an axis in the plane of the circle, lying a distance R from the center of the 
circle. The familiar "doughnut" shape is generated when R > r. Ln this case the large 
circle of radius R lies in the xy plane, and the axis z runs through the center of the 
doughnut "hole". The toric surface used in lenses lies near the outer extremity of such 

10 a doughnut shape, near the point a distance R + r from the center along the x or y axis. 
A segment near this point lying in the xy plane will have a radius of curvature of r + R. 
A segment lying perpendicular to this plane, in the xz plane, for instance, with have 
radius of curvature r. In practice, R is usually significantly smaller than r, so the 
"doughnut" is re-entrant, and very nearly resembles a sphere. The axis along which the 

15 radius is equal to R + r, called the Cylindrical axis (as R + r is called the Cylindrical 
Radius of Curvature) may lie rotated by some angle cp from the Zero axis of reference 
of the SMI. The axis along which the radius measures r (the Spherical Radius of 
Curvature, measured along the Spherical axis) is always perpendicular to the 
Cylindrical axis. 

20 For such a part plots of the coefficients a 0 2, a 2 o, and an as a function of 

position z look somewhat different than in the spherical case, an is still constant, but is 
not equal to zero. The coefficients ao2 and a 2 o will look similar to each other, 
maintaining a fixed separation, but will not, in general, coincide. The two plots will, in 
fact, coincide if cp equals 45°. The values of these three coefficients, along with a 

25 knowledge of the distance from the Surface location at which the coefficients are 
determined, will give the values of r, R, and cp. (The plots of ao2 and a 2 o will still 
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coincide at the Surface location, so there is no ambiguity). The relevant equations are: 



r — 



- + (^20 + O+V( fl 20 - a Olf +Oll) 2 



7? = 



l-W( a 20-O 2 + (^ll) 2 



(ft = — tan" 
2 



<3i 



V^20 



*02 y 



In order to use these equations one need only know the values of the 
coefficients at a single location (whose position relative to the Surface must also be 
known). It is useful to obtain measurements at several positions, however, since then 
one can calculate the functional form via a least-square method and thus compensate 
for random errors in measurement. We find by experience that such random errors can 
be significant, so looking at the trend allows us to obtain a more accurate result than 
we can through measurement at any one position. This is also an advantage of this 
method of performing Shack-Hartmann tests over static methods. This same method 
can be used for the accurate determination of all coefficients up to the order used for 
the fitting of the shift data. Coefficients a nm in which either n or m is odd will be 
constant with position z, while those for which both n and m are even can be expressed 
in the form (constant) - (constant')/z~ (n+ml) . A least square fit to a series of 
measurements of a nm (z) can easily be performed using matrix methods, yielding a 
robust result that is insensitive to noise. 

As another example, the value of the Shape Factor of an aspheric part can be 
obtained by using the values of the fourth order coefficients. The Shape Factor is 
given by k + 1, where k is the "Conic Constant". The value of k + 1 equals 1 for a 
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sphere. If the value is greater than 1 or between 1 and 0 the result is an ellipse. A value 
of k + 1 = 0 gives a parabola, while a value less than zero yields a hyperbola. 

Provided the axis of the part is aligned with the axis of the optical system and 
centered, the Shape Factor k + 1 is given by the fourth order coefficients a^ = ao4 = 
5 a22/2 = a4 by the expression: 

?c+l = 8r 3 a 4 +l 

Here r is, of course, the value of the radius of curvature obtained from crossing-point 
data, or one of the other methods noted above. 

For practical measurements, it may sometimes be necessary to make small 

10 corrections to these mathematically derived formulas. 

Beyond spheres and tones, shapes of increasing complexity can be studied by 
extending the methodology of the invention. By doing this, "atoric" shapes (shapes 
that are toric when observed near the optical axis, but which depart from that shape as 
you get farther from the axis) can be examined as well as more general shapes. There 

15 is a straightforward method for determining the shapes of parts along an arbitrarily 
chosen azimuth, and it is the method used in plotting the radius of curvature maps. 
From a position z along the direction of travel, and the surface parameters a nm can be 
measured. Then the surface can be constructed by adding the sag of a spherical surface 
of radius z (the reference sphere) to the sag determined by these parameters. From this 

20 the radius of curvature can be calculated if desired. 

Thus, the method for measuring spherical parts, and in determining the 
coefficients of toric parts can be extended to more complex parts, measured to higher 
orders. The difference lies in the way the measurements are made. The great strength 
of the invention is the use made of not only the shifts characteristic of the surface 

25 shape at one reference position, but at a series of them. What has been plotted has 

been the values of the second order coefficients a2o, an, and ao2 as a function of z. A 
least-square fit of this data to the theoretical functional form was then performed, and 
results were obtained for the toric parameters that are unusually reproducible. This 
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was possible by using least-square fitting, which allows allows for a sort of averaging 
over positions, "averaging out" the statistical errors at each individual position and 
obtaining a much more stable result. The plots of a 2 o(z), an(z), and ao2(z) might not 
describe precisely defined curves, but the first and last of these can be fit to the 
functional form constant + l/2z (with a different constant in each case) and the middle 
one to constant* and very good fits and results can be obtained. Even use of the 
knowledge that the slopes vs. 1/z must be Vz to perform fine-tuning on the 
magnification of the invention has been used. It is this combination of knowledge of 
the functional form of a nm (z) with least-square fitting that gives us such reliable and 
stable results. 

This procedure can be carried through to surface plotting, radius-mapping, and 
aspheric (and atoric) calculations since there is more information present than use is 
made of in just spherical and toric measurements, and this information resides in 
higher-order terms. By fitting them, via least-square formulations, to their expected 
functional forms, instead of to the variations otherwise seen at a single point, 
advantage of multiple measurements to average out the error can be made. Instead of 
plotting different surface shapes for the results at each position z, a single set of 
parameters is determined which uses is used to plot surface variation, radius map, etc. 
for all positions z. One practical result is that the maps and results do not change as 
when moving through z - there is a single, well-defined shape and set of maps. These 
results are much more reproducible and stable. 

To accomplish this, the functional forms that are to be expected need to be 
known or estimated in a rational way. On the surface, this seems difficult. Surely, it is 
different for a spherical part than for an aspheric part than for a toric, or an atoric part. 

However, the answer is surprisingly simple. The expected forms for all the 
higher order terms are extremely simple, and can be written down without any special 
knowledge of the expected shape of the part. Using the expected forms written below, 
it is a simple matter to perform least-square fitting over a range of z values and obtain 
averaged results for further calculation. 
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The key to this lies in realizing that the shape the surface actually has is a 
constant thing. It is entirely independent of the position z at which the measurement 
takes place. This is true if the surface is spherical, aspheric, toric, or anything else. The 
coefficients a nm will be functions of z, but when you calculate the sag of the surface 
(by adding the sag of the reference sphere of radius z to the sum of a nm x n y m ) J you 
should end up with the same shape, since the surface always has the same shape, 
regardless of where you are standing when you measure it. 

The sag of the reference surface measured at a position z (and presumed to be 
equal to zero at the origin, x = 0, y = 0) can be written in the following way (as we 
have shown often enough): 



Sag(z) = z-^jz 2 -(x 2 +y 2 ) 



Sag(z) = 




c / x x 2 y 2 * 4 

Sag(z) = — + ^ + — j 
2z 2z 8z 



x 2 y 2 y< 



4z J 



8z J 16z 3 



•+ — ^+ — =^-+- 



16z 5 



16z 3 



16z 5 



It will be obvious that the only way that the surface shape obtained by adding 
this sag to the sag made up of the terms a nm x n y m can be independent of z, regardless of 
the shape of the surface, is if the terms a nm can be written in the following forms: 

a 00 (z) = const 
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\*4 



ft i 

ri 



a 1Q (z) = const 
a 01 (z) - const 

1 

fl 9n (z) = 

2z 

a n W = const 

a 02 (z) = const - ^ 
a 30 (z) = con^ 
a 21 (z) = const 

15 a i2( z ) = const 

a 03 (z) = const 

1 



10 



20 



a 40 (z) = const 



8z 3 



<2 31 (z) = const 

1 



a 22 (z) = cows£ - 



4z 3 



25 a 13 (z) = const 



1 

a 04 (z) = c<ms£ 



8z 3 



In the above it should be evident that the constants are generally different. Frequently 
30 most of them are equal to zero. The fact that the coefficients should have the above 
forms receives some verification from the theoretical forms for the coefficients of a 
toric part as given. 
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These coefficients are all of the form listed above. (The third order coefficients 
^30? a 2 i, a 12 , and a 0 3 are all equal to zero, which is a constant). It is important to realize 
that atoric parts are also of this form. One can convince oneself that spherical parts 
are, as well, and even perfectly flat parts. 

This brings up the question of how to determine the coefficients. This should 
be done in a manner similar to the second-order coefficients we determine for toric 
measurements. We perform a least-square fit in which we force the slopes of the 
relevant parameters to be the theoretically fixed values. Most of the coefficients can be 
calculated by simply taking the average over the values obtained at each position z. 
This will be true for all terms that do not show up in the expansion of the reference 
sphere sag. In other words: 



j nm V n J 



l 4- 



average 



This will hold true for a 0 o, ai 0 , an, a 3 o> a 2 i, ai 2 , a 0 3? a 3 i, ai 3 , etc. 

Now, you can determine the second order coefficients from: 
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We then get a 20 (z) = offset 20 + slope/z and a 02 (z) = offseto2 + slope/z. But if we are 
confident of our alignment and magnification factors, then the slope is known to be - 
1/2, and we can determine the offsets from the much simpler formula: 




n 



{offset o 2 ) 



1 



average 



n 



In fact, we can say that if we know that the slope of a term should be -l/bz c , then we 
can write the expression for the average value of the offset as: 



Once these coefficients have all been determined, all that must be done to 
calculate the shape of the surface is to choose a point z and add the sag of the reference 
surface at a point z to the sum of the terms a nm (z)x n y m . Because of the restrictions we 
have placed on these terms, you are guaranteed that the surface thus generated will be 
independent of z, at least to the highest order for which you have generated terms. In 
fact, an alternative way to calculate the surface profile is simply to add up the sum of 
a n mX n y m for the constant terms (a 0 o> aio, a 0 i, an etc.) and the sums of the "offset" parts 
of the varying terms — i.e. offset2ox 2 + offseto 2 y 2 + offset^x 4 + ... This procedure can 
also be used in calculating derived quantities, as well. Calculation of aspheric 
coefficients and shape factors can still be carried out as before, but we now use these 



{offset nm ) 



average 



= a nm( Z k) + 
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"scan-averaged" values, rather than values obtained only at one position z. The result 
should be stabler, more accurate values at all locations. 

The invention can also be used to determine the amount of offset in lathe- 
turned spherical parts. Manufacturers are using lathes to turn spherical parts, but if the 
5 center of the sphere which the computerized lathe tool is moving along does not 

coincide with the axis of the lathe, then there will be a problem. Instead of generating 
a sphere the surface will be something else (a variety of torus, but still very different 
from the shape usually called by that name). An axial slice through such a part will 
either resemble a gothic arch or will have a "cusp" in the center, depending upon 
1 0 whether the offset is positive or negative. 

Such an offset can be determined from the following analysis. If we take z as 
the coordinate along the part (and lathe) axis, and r = (x 2 + y 2 ) 1/2 as the radial 
dimension, then the "Sag" z of a spherical surface of radius R is given by: 



15 



= R-^R 2 -r 2 = R< 
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If the axis is offset by a distance g, then the sag becomes: 



z = R 



1- 



L R J 



The difference in sag between the surface of such a decentered part and a perfectly 
20 centered Reference Sphere of radius Ro is given by Az: 



Az = Rq< 
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If we expand the terms containing the square roots we obtain: 



Az = (R 0 -R)-R 0 
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If we further assume that g « r (as we expect), then we can further simplify this to: 



Az 



^ r 2 (R-R 0 ) + 2grR 0 r A (r 2 - ^ 3 ) + 4gr 3 j R 0 3 r 6 (r 5 -B* ) + 6gr 5 *Q 5 ^ 
2^ + 8 (^) 3 + le^) 5 



If we let g go to zero, then the expression becomes the difference in sag between two 
spheres of different radius (R and Ro), which isn't interesting, and is an old result. On 
the other hand, if we set R = Ro we get the sag between the lathe-turned part and a 
reference sphere of the same radius, which has a very simple and revealing form: 



Az = g 



R 



+ - 
8 



The sag is proportional to the offset g, and is expressed solely in terms of the 
odd powers of (r/R). To lowest order, the error is a conical error gr/R . This is not 
surprising - the tip of that gothic window or that dimple is simply a small cone. This 
is, at least for small values of r, the leading term. As we go further from the axis, of 
course, the higher odd terms will dominate. 
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How do we express that conical term in our present system. We have been 
expressing the departure from the spherical reference in terms of polynomials in x and 
y (or, alternatively, as Zernike polynomials). We express these in the general form: 

But our present expression is in terms of r =(x 2 + y 2 ) 172 , which cannot be expressed as 
a sum of powers of x and y. It is true that if you restrict yourself to values very near 
the x-axis you can approximate r in powers of x and y. (In fact, this is true if you 
restrict yourself to any line with a well-defined relationship between x and y), but that 
is not our case - we consider spot shifts and surface fitting over a large region, 
including a very large area that is nowhere near the x- axis. 

What we want is an approach that uses our current setup and software for a 
quick and simple solution. What if we were to use our present system to map the 
variation in sag from that of a target sphere as a function of radius in terms of the even 
powers we can use (since we can't use odd powers of radius). The simplest way to do 
this is the most direct. It will be described below in detail. 

At any position in the vicinity of the Center of Curvature position (CoC) our 
knowledge of the reference sphere radius (which is equal to the distance we are at from 
the Surface position) enables us to translate the surface deviation h(x,y) into the true 
surface shape by simply adding the reference sphere sag (given by R Re f(l-(l- 
(r/R Re r) 2 ) 1/2 ), where R Re f is the reference sphere radius and r = (x 2 + y 2 ) 1/2 ). We 
simply perform a least-square fit of this true surface shape to that of an off-centerlathed 
piece, using as our variables the offset g and the true radius of curvature of the surface 
RTrue- If the user can supply Rtfuc so much the better - the search for the minimum of 
Chi Squared becomes a one-dimensional search. 

In other words, the software searches for the value of g (and, if necessary, 
RTrue) that minimizes the expression: 
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X 2 = ^{SagiXi , y f , R Ref ) + h(x, yj,R Ref ) - f(x, ,, yj , R True , g)} 

Where: 



Sag(x t , yj ,R Ref ) = R Re / 



h(*i >yj ,RRef ) = X a nm ( R Rc f ) x i 



f( x i >yj>g> R True ) - R True 



The equations are not linear, so the minimization cannot be performed using matrix 
methods. The program must perform a regular search for the minimum of the chi 
squared function, using one of the many methods available. 

This formulation assumes that the part is properly aligned and centered. This 
method can be modified to accommodate decentered and tilted parts, but it would be 
easier if the parts were properly positioned. 

Once the coefficients describing the surface have been obtained using the above 
formulas they can be used to construct an accurate three-dimensional map of the 
surface from which any desired quantity can be derived. This includes, but is not 
limited to, determination of the radius of curvature of a spherical part, determination of 
the shape factor of an aspheric part, determination of the toric parameters of a toric 
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part, determining radius of curvature along selected meridians of a part, determining 
aspheric coefficients along selected meridians of a part, determining the "best fit" 
sphere, parabola, hyperbola, or other defined form to the surface. It can be used to 
determine lateral offset and tilt of a preselected form. It can be used to determine if an 
5 offset exists between a programmed lathe tool and the lathe axis by measuring the 
parameters of the fit. Knowledge of the surface profile as determined by the device 
combined with mathematical models for the expected shapes allow many 
determinations of desired measurements by those skilled in the art. 

The foregoing algorithms for extracting surface and other information from the 

10 measured deviated spot locations can be readily programmed in a variety of program 
languages along with providing the graphical, tabular, and other data values. Fig. 10 
shows, for example, one screen interface that shows a three-dimensional surfaced map 
along with various dialog boxes that provide the user with the ability to select options, 
operations, part description, and display types. 

15 Figs. 1 1A-D and 12A-D show contour and three-dimensional maps of a surface 

as a function of defocus position, Figs. 13 and 14 are output screens illustrating 
coefficient plots versus position of the wavefront origin with respect to a test surface 
for a toric and sphere, respectively, and Fig. 1 5 shows one possible format for 
displaying toric and/or spherical information. 

20 One possible general flow chart consistent with the inventive method is given 

by the following: 

Input the type of measurement desired: 

a. ) Generate Surface Profile only 

b. ) Measure Radius of Curvature of Spherical Part 

25 c.) Measure Spherical Radius of Curvature, Cylindrical Radius 

of Curvature, and Angle for Toric Part 

d. ) Measure Shape Factor or Conic Constant of Aspheric Part 

e. ) Measure Aspheric Coefficients for Aspheric Part 

f. ) Other 
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1. ) Input Target Radius, if known, or indicate that Radius is unknown. Input 

other data (toric angle, etc.) if relevant. 

2. ) Perform measurement of spot deviation in the vicinity of the Surface 

location (i.e., with the translation stage positioned so that the focal point of 
the objective lies near the surface of the part under test. Locate the position 
at which the spot locations differ from the reference positions by at most a 
linear shift. Designate this the "zero" position for the stage, and refer all 
other positions to this one. If no "surface" location can be found produce 
an "error" message. This indicates that the surface finish of the part is too 
poor for measurement, or that the surface at the center is tilted too far from 
normal to the optical axis, or that no part is present. 

3. ) Go to the vicinity about one radius of curvature from this location. If the 

radius is not known, search for the location at which one again finds a 
disposition of spots similar to that found at the reference location or at the 
Surface position. If no such position can be found, report this in an "error" 
message. It indicates that the radius of curvature of the part exceeds the 
range of the machine, or else the finish of the part is too poor to allow 
measurement or there is no part or there is extreme misalignment. 

4. ) Perform at least one measurement but preferably a series of measurements 

at the "Center of Curvature" (CoC) location, similar to that performed near 
the "Surface" location. Measurements are preferably performed at several 
locations in the vicinity of the CoC. At each position the shifts in the spot 
locations are measured. From these shifts and the known distance z from 
the surface (and the known parameters of the SMI device) the coefficients 
aio, a 0 i, etc. of the polynomial describing the deviation of the surface from 
a reference sphere of radius z are determined. 

5. ) Surface maps of the surface deviation at selected positions are generated 

and displayed. 

6. ) If a specific Measurement has been requested the software uses the relevant 

formulae to calculate the needed quantities from the positions z and the 
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coefficients a nm (z). For spherical parts the Radius of Curvature can be 
calculated. For Toric parts the two Radii of Curvature and the angles the 
axes of the radii make with the reference axis of the SMI are calculated. For 
Conic Section parts the Conic Constant or Shape Factor is calculated (as 
well as the base Radius of Curvature). Other measurements may be added 
to this list. Least square algorithms can be applied to fit predetermined 
analytical functions, spheres, aspheres, torics, or the like, to the measured 
surface or wavefront shape. 
Those skilled in the art will be able to imagine other variants based on the 
teaching of the disclosure. Some of those variants are shown in Figs. 16 through 21. 

Fig. 16 is a diagrammatic elevational view similar to Fig. 5 with the addition of 
a stationary positive lens 82 located forward of the microscope objective 70 to enable 
measurement of long radius parts where parts from Fig. 5 that reappear here in Fig. 16 
bear the same numerical identification they had in Fig. 5. If lens 82 has a focal length, 
f 2 , then the radius of a part as indicated at 84 is given by: 

R=f 2 2 /x 

where x is the displacement of the microscope objective 70 required for collimation 
off of the reflective part 84. If x is positive, then R is convex, and if x is negative, 
then R is concave. 

Reference is now made to Fig. 17 which shows a diagrammatic elevational 
view of an embodiment of the invention that is similar to Fig. 5 except for the absence 
of the objective which has been removed to enable measurement of flat surfaces that 
are at least partially reflective. With the objective removed, a flat object such as that 
designated at 86 can be tested for the flatness of its partially reflective surface by 
causing it to be illuminated with a plane wavefront 75 that now becomes the reference 
against which any wavefront distortion caused by the test surface may be compared as 
before. Again, repeated parts are identified here with their previous number. 
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Fig. 18 is a diagrammatic elevational view similar to Fig. 17 except for the 
presence of an additional mirror and a relay system for measuring the wavefront error 
of transmissive filters over their bandpass. As can be seen here, a relay system 
comprising lenses 88 and 90 is placed in the intervening space between beamsplitter 80 
5 and a bandpass filter 92. Plane reference wavefront 75 passes through bandpass filter 
92 after which it impinges on flat mirror 94 from which it reflects to travel back 
through bandpass filter 92, now having made a double pass through it, enters the relay 
system, and emerges as a more or less distorted wavefront which enters the camera 76. 
Because the filter transmits only the light within its bandpass, this embodiment 

10 measures the wavefront error of the filter only within that bandpass. Thus, the use of a 
broadband source instead of a spectrally narrow source has the advantage of being able 
to test filters whose bandpasses can span a broad spectral range. 

Fig. 19 is a diagrammatic elevational view of an embodiment of the invention 
similar to that shown in Fig. 1 7 with the addition of a beam expansion section to 

15 measure large flat surfaces. The beam expansion section comprises two positive lenses 
96 and 98 that operate in conjunction with one another to enlarge the diameter of the 
plane reference wavefront so that it can cover a larger area of a reflective part to be 
measured as, for example, the flat designated at 99. Alternative telescope geometries 
can also be used. 

20 Fig. 20 is a diagrammatic elevational view of an embodiment of the invention 

similar to that of Fig. 18 except for the absence of a relay section and, like the 
embodiment of Fig. 18 measures the wavefront error of transmissive filters over their 
bandpass. This embodiment is suitable for use where a relay section is not needed to 
image the part under test on to the lenslet array. 

25 Reference is now made to Fig. 21 which is a diagrammatic elevational view of 

an alternative embodiment of the invention that uses an intervening telescope between 
the objective and CCD array to make full use of the available numerical aperture (NA) 
of the objective. Here, there is an objective 100, a CCD array 102, and crossed 
lenticular lens 104, all of which operate as before. The optical path between objective 

30 100 and CCD 102 includes an intervening telescope comprising lenses 110 and 112 
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and a pair of mirrors 106 and 108 for folding it for vertical compactness. The purpose 
of the telescopic lenses 110 and 1 12 is to optically scale the size of the CCD and 
lenticular array so that they fill the pupil of objective 100 to make full use of its 
numerical aperture (NA). As a convenience, a light source 1 16 is introduced in the 
5 telescope via beamsplitter 114 that is conjugate to a spatial filter 120 that operates to 
eliminate stray light. 

Thus, it will be appreciated that spot sizes, CCD resolution, pixel saturation 
levels, objective NA and magnification, and pinhole size are among design parameters 
that may be varied. In addition, it will be appreciated that scaling optics may be placed 

10 between the microscope objective and the lenticular screen or the image plane of the 
lenticular array and the CCD to reimage the array of spots onto the CCD to rescale the 
relative sizes of the spots with respect to the pixels. Trigonometric polynomials or 
other series approximations for surface shape may be used. It will also be appreciated 
that the thickness and power of an element may be derived by examining both front 

15 and rear surfaces. In addition, the use of a strobe light or pulsed source may be used in 
place of a continuous source to study the dynamic properties of surfaces and 
transmitting components and the apparatus and test samples may be moved relative to 
one another to scan the properties of the surface in, for example, a web or scanning 
process. It will also be clear that the two-dimensional CCD array may be replaced by a 

20 one-dimensional array that scans across the focal plane of the crossed lenticular screen 
array or its equivalent. Accordingly, it is intended that these and other modifications 
and variations be within the scope of the invention. 
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